Small Tj coherence peak near T c in unconventional BCS superconductors 
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It is usually believed that a coherence peak just below T c in the nuclear spin lattice relaxation rate ly 1 in 
superconducting materials is a signature of conventional s-wave pairing. In this paper we demonstrate that 
any unconventional superconductor obeying BCS pure-case weak-coupling theory should show a small Tj -1 
coherence peak near T f , generally with a height between 3 and 15 percent greater than the normal state Ty 1 at 
T c . It is largely due to impurity and magnetic effects that this peak has not commonly been observed. 



PACS numbers: 

It is well known 1 1] that conventional superconductors obeying 
BCS weak-coupling theory generally show a large coherence peak in 
the nuclear spin lattice relaxation rate Ty 1 below T c . This is a direct 
consequence of the large quasi-particle density-of-states found for 
E > 1 .0 at the gap edge. It is commonly believed, and experimentally 
largely true, that unconventional superconductors show no coherence 
peak near T c . However, a simple quantitative argument shows that if 
the BCS weak-coupling pure case theory applies, a small coherence 
peak below T c must exist. 

Consider the BCS weak-coupling equation for the nuclear spin 
lattice relaxation rate T^ 1 for unconventional superconductors: 

{TyT)- l /{T x T)- T x =Tc = ^dEN 2 (E)sech 2 (E/(2T))/2T (1) 

The sech 2 (£/(2r)) acts as an attenuation factor and dominates the 
low-temperature Tj , yielding exponentially activated behavior for 
s-wave superconductivity and power-law behavior for unconven- 
tional superconductivity. However, as T — > T c all of the structure 
in N(E) (i.e. DOS different from unity) is shifted to lower energies, 
since N(E) = N(E/A) and A — > 0. See Figure 1 for a depiction of this 
behavior for the 3-d 3 -He A-phase order parameter A(k) = A sin 9. 
All of the structure in N(E) falls in the region where E/2T C 1, so 




FIG. 1: Diagram depicting attenuation factor and quasi-particle 
density-of-states for T ~ 0.98r c ~ 0.5A ,A = 0.2A , 3-d 3 -He A- 
phase order parameter. 

that for this structure the exponential damping factor becomes essen- 
tially irrelevant. Define F(E) = N(E) — 1. Note that by the density- 



of-states sum rule, J ( ™ F(E)dE = 0. Now substituting in for N(E) we 
find 

f dEn+2F(E)+F 2 (E))sech 2 (E/(2T))/2T (2) 
Jo 

The first term trivially yields 1. The second term, 

/ °° dE 2F (E) sech 2 (E/(2T))/2T can be evaluated by noticing 

that as A — » 0, F(E) is only significantly different from in regions 

where the argument of the sech 2 is small, so that to an excellent 

approximation near T r this integral is equal to J{fE(E)dE = 0. 

The final term, dE F 2 (E)sech 2 (E/(2T))/2T, is positive, and 

so in the immediate neig hborhood of T c , (7i T)~ l /(T\ T)T^_ T ^ > 1, 

implying the existence of a peak. It is this redistribution of N(E) 

away from an energy-constant (= A^o) density-of-states, represented 

by F 2 (E), that is responsible for the peak in Tf near T c . The larger 

this effect, the larger the peak. 

This redistribution is intimately tied in with the nodal structure of 

A(k). This can be seen directly from the BCS expression for the 

density-of-states N(E/A) = N(x) = Re( ■'" = ), where (...) de- 

v f~ 

notes an average over the Fermi surface and f contains the angular 
dependence of the order parameter, (i.e., A(k) = Ao/(k)). The con- 
tribution of the nodes is most easily parametrized by < f 2 >, with 
larger values indicating less nodal order parameters. For example, 
an s-wave order parameter has < f 2 > = 1 , while a 2d d-wave order 
parameter (containing line nodes) has < f 2 > =0.5. Gap functions 
/ with larger < f 2 >, indicating effectively small or absent nodes, 
have a comparatively smaller region of phase space contributing to 
the integral, for x < 1 . These gap functions will therefore show de- 
pleted low-energy density-of-states, and by the sum rule must have 
enhanced spectral weight in the peak at E = A. Both effects will tend 
to enhance the T^ 1 peak just below T c . 

These behaviors are illustrated in Figure 2, which depicts 
densities-of-states and Tf for a series of 3-d order parameters 
A(k) = 1 — cos" (9), with 9 the polar angle. As n increases, the low- 
energy DOS is depleted and the coherence-peak DOS enhanced, with 
a concomitant increase in the Tf peak near T c . For these cases, 
< f 2 > increases monotonically from ~ 0.533 for n=2 to 0.866 
forn=10. 

Results. - In Figure 3 are depicted the coherence peaks near T c for 
several unconventional order parameters: d v 2_ v 2-wave, the 3 -He A- 
phase order parameter (for which A(k) = Ao sin 9), and the p-wave 3- 
dimensional order parameter A(k) = Ay cos 9, as well as the quasipar- 
ticle density-of-states for these order parameters. Note that a small 
peak just below T c is evident even for the last order parameter, whose 
quasiparticle density-of-states shows no divergence at E = A. The 
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the elliptical function: 



FIG. 2: The densities-of-states and Tj 1 for several 3D order param- 
eters A(k) = 1 — cos"(0) are shown. 
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Now, to work out an analytic form for the peak in Tj just be- 
low T c , we must compute /J° dEF 2 (E)sech 2 (E /(2T))/2T, where 
F(E) = N(E) — 1. For T sufficiently near T e , F(E) only varies from 
zero in a region where sech 2 (£'/2r) is essentially unity, as described 
at the beginning of this paper. Making the substitution E = Ax, and 
taking sech 2 (£/2r) as 1, we find for the above integrals, where yis 
the Catalan constant= 0.915 . . . 



3D line node 




FIG. 3: The predicted nuclear spin lattice relaxation rates (TiT) 
near T c and quasiparticle densities-of-states are shown for three order 
parameters: 2D d-wave, the A phase of 3 //e, and a 3D line-node 
model. 



analysis of the preceding paragraph demonstrates that it is largely the 
second moment of the DOS around an energy-constant DOS (~ No) 
that produces the small peak just below T c . While a large or diver- 
gent DOS at E — A clearly enhances the coherence peak near T e , it 
is not necessary for the formation of a peak. 

For the three cases described above, it is possible to derive an 
analytic expression for T^ 1 just below T e and compare with the nu- 
merical results. The quasiparticle density-of-states for each of the 
three order parameters can be computed analytically, and one finds 
the following well known results QSHl' where x = E/A and K is 



d-wave: ~ 0.4743A/2F (8) 
, it 2 

3 -HeA-phase : — A/2r ~ 0.822A/2r (9) 

,2 1., 

3D line node: (-y- -)A/2T ~ 0.2773A/2r (10) 

In other words, very near T c (T > 0.995r f ) we can express the ra- 
tio (T] r)~' /(T\T)p_ T as simply l+aA(T)/T, where a is an order 
parameter-dependent constant, and this expression yields reasonably 
good agreement with the numerical results. In order to better model 
the behavior near T c we have calculated analytically the next order 
term and found excellent agreement, as indicated in the plot below. 
Below 0.98 T c this approximation becomes less accurate. 




FIG. 4: The analytic and numerically computed Tj~ for the three 
order parameters previously discussed are shown. Symbols indicate 
the analytic result and the lines the numerical result. 

Discussion.- The foregoing analysis shows that BCS weak- 
coupling pure-case unconventional superconductors should exhibit 
a small T^ 1 coheren ce peak j ust below T c . Yet a literature su rvey on 

^topiEJi^ b e m m s m nam 

I22L I23L 124. 125I l2(i I27II has uncovered just two materials - CePt3Si 
and (TMTSF^PFg-which show such a peak. The question therefore 
arises as to why such peaks are not commonly observed. 
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To address this question, we have conducted an analysis of the 
effect of resonant impurity scattering upon this T^ 1 peak, for a two 
dimensional d-wave order parameter. It is well known that such im- 
purity scattering truncates the DOS peak at E = A and in addition 
can generate substantial low-energy density-of-states. Both of these 
effects would tend to reduce the size of the coherence peak. It turns 
out that for these reasons the appearance of this peak is extraordi- 
narily sensitive to impurity scattering. Depicted below are four Tj~ 1 
curves for d-wave superconductivity: zero impurity scattering, and 
three cases of small impurity scattering: r/Ano = 0.01, 0.02 and 
0.03. Within the unitary limit this last concentration is roughly 7 
percent of the critical impurity concentration required to destroy su- 



fermion and cuprate materials. Indeed, the organic superconductor 
(TMTSF^PFs has already shown a small peak below T c l27ll . 

To summarize, here we have demonstrated that any unconven- 
tional superconductor obeying BCS pure-case theory should show 
a small coherence peak in the nuclear spin lattice relaxation rate Tj" 1 
just below T c . It is likely due to magnetic and impurity effects that 
this peak has not generally been observed in unconventional super- 
conductors. 
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FIG. 5: The numerically computed Tj for no impurity scattering, 
and for r/A 00 = 0.01, 0.02 and 0.03 are shown. 

perconductivity 1 28], and would result in a depression of T c of this 
order. As is clear from the plot, the height of the peak is greatly di- 
minished even by the rather low impurity scattering 1 3] rates modeled 
here. 

Given that materials in the unitary limit typically have supercon- 
ductivity destroyed by an impurity concentration on the order of a 
few percent ] 29] , the foregoing analysis indicates that an impurity 
concentration of just 0.25 percent is sufficient to largely destroy this 
peak. Such a concentration is well within the range of observation 

In order to observe this peak samples of the highest possible qual- 
ity are clearly essential, with impurity concentration less than 0.1 
percent. It would also be advantageous to perform low-temperature 
specific heat measurements on the same samples as this would al- 
low accurate assessment of the prediction of a finite relaxation rate 
at T=0, via a measurement of the residual density of states. 

An additional effect complicating the observance of this peak 
is the frequent occurrence of magnetism in the heavy-fermion and 
high-T c cuprate materials upon which most of the measurements 
have been performed. The combination of magnetism and the ex- 
treme sensitivity of this peak to the presence of impurities in the 
unitary limit make its observation difficult in the heavy-fermion and 
cuprate superconductors. However, the recently discovered non- 
centrosymmetric superconductor Li^PhB |3lll shows no signs of 
magnetism or strong electron correlation [32], and appears to be un- 
conventional on at least one band, based upon magnetic penetration 
depth data 1331 . This material may therefore be an ideal material in 
which to search for this small T^ 1 peak. Another possibility for ex- 
periment is the class of organic superconductors, which may not nec- 
essarily have the sensitivity to impurities characteristic of the heavy- 
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